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Abstract. We analyze in some detail the recently discovered velocity quantization 
phenomena in the classical motion of an idealized one-dimensional solid lubricant, 
consisting of a harmonic chain interposed between two periodic sliders. The ratio 
w = v cm /v ex t of the chain center-of-mass velocity to the externally imposed relative 
velocity of the sliders is pinned to exact "plateau" values for wide ranges of parameters, 
such as sliders corrugation amplitudes, external velocity, chain stiffness and dissipation, 
and is strictly determined by the commcnsurability ratios alone. The phenomenon is 
caused by one slider rigidly dragging the density solitons (kinks/antikinks) that the 
chain forms with the other slider. Possible consequences of these results for some real 
systems are discussed. 
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Figure 1. (a) A sketch of the model with the two periodic sliders (periods a+ and a_) 
and the lubricant harmonic chain of rest nearest-neighbor length <Xq. (b) The finite-size 
phase error A of Eq. |J2) (only the N values where A is smaller than for any other 
N 1 < N are shown) decreases as a function of size N. Note that the error drops fast 
as N~ 1 for golden mean and other quadratic incommensurabilities (r + ,r_ = r?_), but 
slowly as N~ x / 2 for spiral mean and other non-quadratic incommensurabilities. 

energy a + F + 

velocity v a 1 / 2 F]/ 2 mT x l 2 

time ' afF + - 1/2 m^ 

spring constant K a^ 1 F + 

viscous friction 7 a^ 1 ^ 2 F]/ 2 m 1//2 



Table 1. Natural units for several physical quantities in a system where length, force 
and mass are measured in units of a + , F+, and m respectively. 

1. Introduction to the model 

The sliding friction of hard incommensurate crystals has stimulated for over a decade the 
study of non-linear classical dynamical systems - belonging to the class of generalized 
Frenkel-Kontorova (FK) models jl] - which can show a rich variety of dynamical 
behavior, ranging from regular to chaotic motion [2J. In particular, in presence 
of a uniform external driving, the length-scale competition between different spatial 
periodicities can give rise to intriguing frequency incommensurations and phase locking 
phenomena. In this paper we discuss and analyze in depth the striking exact velocity 
quantization phenomena recently reported in a one- dimensional (ID) non-linear model 
inspired by the tribological problem of two sliding surfaces with a thin solid lubricant 
layer in between [3J 0] . The model "layer" consists here of a chain of N harmonically 
interacting particles interposed between two rigid generally (but not necessarily) 
incommensurate sinusoidal substrates representing the two "sliding crystals" , sketched 
in Fig. [Tfe,, externally driven at a constant relative velocity t> ex t- The equation of motion 
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of the z-th particle is: 



mxi = — — [F + sin k+(xi — v+t) + F_ sin k-{xi — v-t)] 
+ K(x i+l + ar<_i - 2xi) - T^X^ ~ u ±) > 



(1) 



± 



where m is the mass of the N particles, K is the chain spring constant, and k± = (2ir) / a± 
are the wave-vector periodicities of potentials representing the two sliders, moving at 
velocities v±. We choose the reference frame where v + = 0, so that t> ext =V- — v+ = u_. 
7 is a phenomenological parameter substituting for various sources of dissipation. 
Dissipation is required to achieve a stationary state, but has otherwise no major role in 
the following. F± are the amplitudes of the forces representing the sinusoidal corrugation 
of the two sliders (we will commonly use F-/F + = 1 but we did check that our results 
are more general). We take a + — 1, m — 1, and F + = 1 as our basic units, so that the 
natural units for all quantities used to characterize the model are listed in Table HJ 

The relevant length ratios [51 E] are therefore r± = a±/ao; we take, without loss 
of generality, r_ > r + . The inter-particle equilibrium length a , not entering explicitly 
the equation of motion (CQ), appears via the boundary conditions, which are taken to be 
periodic (PBC), x N+1 — x\ + N a , to enforce a fixed-density condition for the chain 
PQ, with a coverage r + of chain atoms on the "denser" substrate of periodicity a + . 
Technically, the PBCs can be realized in two different manners. The first (i) is quite 
standard [T], and consists in approximating r + and r_ with suitable rational numbers 
N/N + and N/N- (with N± = [N/r±] meaning the nearest-integer of N/r±). The second 
(ii) uses, instead, a finite iV and machine-precision values of r + and r_. For irrational 
r + and/or r_, both methods introduce a finite-size error 



in the relative phase of the three periodic objects involved in the model. In method (i) 
this error is distributed uniformly through the chain, while in method (ii) it remains 
concentrated where the boundary condition applies, at sites 1 and N. For a generic 
N, A is of order unity, but, as illustrated in Fig. [Tb (reporting the N values where A 
is smaller than for any other N' < N), it can be made arbitrarily small by choosing 
suitably large integers N,N±. In both methods we obtain, as long as A < 0.1, fairly 
stable and consistent results, which converge to the same limit for large N, N±, where A 
vanishes and full incommensurability is restored. In a (r + , r_ —rV) scheme, the speed of 
convergence depends crucially on the nature of the incommensurability of r+. Quadratic 
irrational values of r + are, thanks to their periodic continued-fraction representation [7], 
especially advantageous and yield a rapidly vanishing phase error A oc iV -1 (the example 
quadratic irrationals of Fig. [lb include the golden mean <p = (1 + V5)/2 ~ 1.618, with 
continued-fraction coefficients (1,1,1,...), and (2 + V2)/3 ~ 1.138, with continued- 
fraction coefficients (1, 7, 4, 8, 4, 8, 4, 8, • • •)). By contrast, for non-quadratic irrationals, 
and/or for general r_ ^ r 2 + , the phase error vanishes much more slowly A oc N~ 1 / 2 (the 
examples of Fig. [lb include the spiral mean o ~ (1, 3, 12, 1, 1, 3, 2, 3, 2, 4, ■ ■ ■) ~ 1.325 



A = A, + A 



with A± = 2vr N± 



(2) 
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solution of a 3 = a + 1, and tt/2 ~ (1, 1, 1, 3, 31, 1, 145, 1, 4, 2, • • •) ~ 1.571). Finally, for 
the class (r + , r_ = 6r + /[r + — 1]), with rational 9, the error decays rapidly A oc N~ x , 
regardless of the irrational nature of r+. 

Finite chains with open boundary conditions (OBC) were also considered, with 
similar results briefly discussed below, and dealt with in greater detail in Ref. [8]. 
Similar models were considered in the past [91 [10l [11] missing however the ingredient of 
incommensurability, which is important for the novel features shown here. A previous 
study of the present model [12] achieved the sliding through application of a driving 
to one of the two substrates, via an additional spring. That procedure obscures 
the quantization phenomena which are instead uncovered when sliding occurs with a 
constant velocity v ext . 

Upon sliding the substrates, t> ex t 7^ 0, the lubricant chain slides too. However, it 
generally does so in a nontrivial manner: the time-averaged relative chain velocity w = 
v cm /v cxt , is generally asymmetric, namely different from 1/2. Even more surprisingly, w 
is exactly quantized, for large parameter intervals, to plateau values that depend solely 
on the chosen commensurability ratios. The asymmetrical w-plateaus are generally 
very stable, and insensitive to many details of the model, indicating an intrinsically 
topological nature. We show that they are the manifestation of a certain density of 
topological solitons (kinks) in the lubricant versus the "denser" substrate potential with 
periodicity a+. These solitons are set into motion with the external driving velocity v cxt 
by the sliding "sparser" substrate with periodicity a_. 

2. The plateau dynamics 

We now turn to illustrating the results. We integrate the equations of motion (TjQ) starting 
from fully relaxed springs (x{ = ia , ii = v ext /2), by a standard fourth-order Runge- 
Kutta method. After an initial transient, the system reaches its dynamical stationary 
state, at least so long as 7 is not exactly zero. Figure [2^ shows the resulting time- 
averaged center-of-mass (CM) velocity function of the chain stiffness K for 
four representative (r + ,r_) values: three with r + > 1 and r_ = r\, and one with 
r + < 1. We find that w = v cra /v ext is generally a complicated function of K, with flat 
plateaus and regimes of continuous evolution, in a way which is qualitatively similar for 
different cases. The main surprise is that all plateaus show perfectly flat w values, that 
are exactly constant (quantized) to all figures of numerical accuracy, the precise value 
strikingly independent of K, of 7, of t> ext , and even of 

To explore the origin of the w plateaus, we analyze the dynamics for a large number 
of values of (r + ,r_), and observe that: (i) at least one velocity plateau as a function 
of K occurs for a wide range of (r + ,r_); (ii) additional narrower secondary plateaus 
often arise for stiffer chain (larger K, see Fig. Ek); (iii) the velocity ratio w of the main 
plateau - the first plateau found for increasing K - satisfies 

w = 1 - r- 1 , (3) 

for a large range of (r + ,r_), see Fig. [2b. 
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Figure 2. (a) Average drift velocity ratio w = v cm /v cx t of the chain as a function 
of its spring stiffness K for different length ratios (r+,r_): commensurate (3/2,9/4), 
golden mean (GM) (</>, 2 ), spiral mean (a, a 2 ), r + close to unity (yjj§|j, i^VlOl) ~ 
(1.031,10.36) - dilute kinks, and r + smaller than unity (</> _1 ,</>) - antikinks. When 
relevant, plateaus are labeled by the rational 9 of Eq. ((6|). All simulations are carried 
out in PBC with 7 = 0.1, v cxt = 0.1. The (0, 4> 2 ) plateau value is w = 0.381966 . . ., 
identical to 1 — to eight decimal places, (b) The main plateau speed w as a function 
of r+. 
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interparticle distance 

Figure 3. Snapshots of the distance between neighbor lubricant particles in the chain 
Xi — Xi-i at three successive time frames. All simulations are carried out with 7 = 0.1, 
Woxt = 0.1, r_ ~ 10.36, and K — 10 (inside the main plateau), (a) r + = 1.031 (kink 
density S/a + = 0.031/a+); (b) r+ = 0.995 (anti-kink density \5\/a+ = 0.005/a+). 



Solitons and exact velocity quantization of incommensurate sliders 



6 



As sketched in our recent note [3], we can understand these results as follows. 
Consider for illustration a situation of quasi-commensuration of the chain to the a + 
substrate: r + = 1 + 5, with S <C 1. The slight commensurability mismatch induces a 
density 

Psol = — = (4) 

of solitons (or kinks), essentially substrate minima holding two particles, rather than 
one pp. Assume that the second, less oscillating a_ slider, which moves at velocity t> ext , 
will drag the kinks along: v so i = v ext . If po = l/ a o = r + /a + is the linear density of 
lubricant particles, mass transport will obey v cm p = v so i p so \. This yields precisely 

w _ Vcm _ Psol_ _ _ j 1_ 

v cx t Po r+ r + ' 

Thus the quantized velocity plateaus appear to arise because the smoother slider (whose 
exact period a_ is irrelevant, but is the one further away from the chain's periodicity, 
r_ > r + ) drags the lattice of kinks, of fixed and given density, at its own full external 
speed f ext , as illustrated Fig. [3k. As shown in Fig. [2b, the physics of kink dragging 
is not limited to the quasi commensurate case i5 C 1, but extends to all values of 
\8\ ~ 1, where individual kinks can hardly be singled out. Amusingly, this mechanism 
also works for 5 < i.e. r+ < 1 . Here kinks are replaced by anti-kinks, and we find 
that the chain moves in the opposite direction to the external driving t> ext (w < 0, see 
Fig. [2k for r + = _1 ). Exactly as holes in a semiconductor, anti-kinks (carrying a 
negative "charge" ) dragged at velocity +v ex t effectively produce a backward net motion 
of the lubricant layer. This "upstream" motion of the anti-kinks (regions of increased 
inter-particle separation, and of decreased local particle density) is illustrated in Fig. [3b. 

Before moving on, we must explain what makes the chain select the substrate with 
which kinks are formed, as opposed to the other substrate merely acting as a moving 
brush dragging the kinks along. As it turns out, kinks are formed relative to the nearer- 
in-register substrate, which by construction has here the period a + . The density (j^) of 
kinks is such that they can be either commensurate or incommensurate relative to the 
periodicity a_ of the second substrate. The soliton lattice is commensurate to the a_ 
substrate whenever the ratio of a_ to the average inter-kink distance is a rational 
number 9, i.e. when 

' A 0. (6) 



r + — 1 

We find that this condition, for 9 = 1, is particularly beneficial to formation of a large 
and stable main plateau. Indeed (r+, r_) = (0, x^) = (0, 2 ) does satisfy Eq. [6] for 
9 = 1, and its i^-plateau covers almost three decades, as shown in Fig. [2k.. We find 
similarly wide plateaus for (r + , ;rzj), with r + = a, r + = n/2 and r + = 7T 1 / 2 . 

2.1. Individual particle motion and CM periodic dynamics 

The motion of the individual particles is also very instructive. Figure H] shows the 
time evolution of the velocity of a single chain particle Xj, and of v cm , for a value of 
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Figure 4. Time evolution of a particle velocity Xi, and of the chain CM velocity 
v cm (fluctuations rescaled by a factor 50), for the GM (a) and SM (b) cases of Fig. [2] 
Amplitudes of the Fourier spectrum of ii(t) (c) and of v cm (t) (d) for r = a. Individual 
particle spectra have identical amplitudes, and differ only in the phases, which leads 
to a remarkable cancellation in the w cm power spectrum. Here K = 1, 7 = Q.l, and 
v cxt =0.1. 



K inside a plateau, K — 1, for r + = <p and r + = cr of Fig. [2j Here two clear kinds 
of behavior emerge. The single-particle motion of each particle in the chain is, during 
motion in the GM plateau, perfectly time-periodic. A similar periodic dynamics is 
found, in appropriate regimes, for all rational and quadratic irrational (r+, r_ = ri) 
values which we tested. Similarly, periodic motion emerges on the plateaus of (r + ,r_) 
satisfying Eq. Q, for arbitrary r + . On the contrary, single-particle motion in the (cr, a 2 ) 
spiral mean (SM) case is definitely not periodic. The Fourier spectrum of the particle 
motion, shown in Fig. Ht, confirms that SM sliding yields only quasi-periodic orbits with 
two prominent incommensurate frequencies /+ and /_, which can be interpreted as the 
average frequency of encounter of a generic particle with the two substrate periodic 
corrugations, /+ = v cm /a + and /_ = (v cxt — t> cm )/a_ respectively. There is however a 
complete and exact phase cancellation between the Fourier spectra of different chain 
particles (all having the same amplitude spectrum, with different phases), giving rise to 
a strictly periodic motion of the chain CM even in the SM case. Periodic CM velocity 
oscillations around an exactly quantized drift velocity is a common and generic feature of 
all plateaus in the chain dynamics. These periodic oscillations can be understood as the 
solitons moving at velocity t> ext encountering the periodic Peierls-Nabarro potential pQ of 
period a+. The corresponding frequency of encounter t> ext /a + is clearly visible in Fig. HU. 
We observe that the Peierls-Nabarro barrier vanishes (yielding a strictly constant v cm ) 
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for all cases where the particle motion is periodic, particularly all (r + ,r_) satisfying 
Eq. ([6]). When such a periodic single-particle motion occurs for a a_-commensurate 
soliton lattice, Eq. (jBJ), the frequency ratio 



As this condition is equivalent to Eq. (jHJ), this analysis attributes to the rational 9 the 
double significance of (i) the coverage fraction of solitons on the substrate of periodicity 
a_, and (ii) the ratio between the average frequencies /+ and /_ of encounter of the 
features of the two substrates. 

Low driving velocities t> ex t are beneficial to the appearance and to the width of 
velocity plateaus. For increasing t> ext , the plateaus shrink and eventually disappear, 
still remaining exact while they do so. The critical t> ext where the plateaus eventually 
disappear depends on K, but is usually smaller than unity, for the parameters of Fig. [2j 
Fully commensurate 9 = 1 plateaus such as that of the (0, (p 2 ) case are especially wide 
and robust against an increase of v ext (for the parameters of Fig. [2] and K = 4, up 
to v cxt ~ 1.5) and other perturbations. Within the class of (r + , rj_) length ratios, 
r + = cf) appears therefore, in the present dynamical context, as the "most commensurate" 
irrational, at variance with static pinning in the standard FK model, where the opposite 
occurs [TJ [T3] . 

All the above is for PBC. As anticipated, OBC simulations [8] show however that the 
PBC used are not crucial to the plateau quantization, which occurs even for a lubricant 
"patch" of finite and not particularly large size N, perhaps such as a hydrocarbon 
chain molecule would be in ID, or an interposed graphite flake in two-dimensions (2D). 
The finite lubricant size does not remove the plateaus. On the contrary, it permits and 
favors the development of robust velocity plateaus as a function of stiffness. The plateau 
values are not identical to those of the infinite chain with PBCs. The reason for that 
is that finite size allows for an overall chain- length re- adjustment which spontaneously 
promotes single-particle periodic oscillations, effectively satisfying Eq. (jSJ). 

2.2. Dynamical incompressibility and tribological dissipation 

Returning to the general infinite chain case, the finding of exact plateaus implies a 
kind of "dynamical incompressibility" , namely identically null response to perturbations 
or fluctuations trying to deflect the CM velocity away from its quantized value. In 
order to probe the robustness of the plateau attractors, we can introduce an additional 
constant force F acting on all particles in the chain, whose action is to try to alter 
the force-free sliding velocity away from its plateau value. As expected, as long as F 
remains sufficiently small, it does perturb the single-particle motions but has no effect 
whatsoever on w, which remains exactly pinned to the attractor (t> cm = i> p i a teau)- The 
plateau dynamics is only abandoned above a critical force F c . This dynamical depinning 
takes place through a series of type-I intermittencies [14], as shown in Fig. [5b where 
v cm (t) is plotted against a slow adiabatic ramping of F. 




9. 



(7) 
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Figure 5. Dynamical dcpinning force F c for the GM plateau (a), extracted by a slow 
adiabatic increase of the force F applied to all particles, until intcrmittencies appear, 
signaling collective slips (b) . The critical behavior of the average velocity for K — > K c 
(c,d), showing the typical square-root singularity (dotted line) associated to type-I 
intcrmittencies. Here 7 = 0.1, and i> ex t = 0.1 were used. 



A precise value of F c can be obtained by ramping F with time with a gentle enough 
rate of increase or, alternatively, by a Floquet-Lyapunov linear stability analysis. This 
procedure is standard [15] and amounts to studying the eigenvalues of a 2N x 2N matrix 
R whose columns are the linearized time-evolutions, over one period T, of the standard 
basis vectors e, = 5ij in phase-space. The motion is stable if the (complex) eigenvalues 
Vi of R are smaller than unity, \v<\ < 1; intermittencies of type I arise when the largest 
(in modulus) eigenvalue of R, ^ max , reaches the point 1 in the complex plane, z/ max — > 1, 
which is indeed what we find at the edge of the GM plateau. The value of F c is a 
function of the parameters, and F c vanishes linearly when K approaches the border K c 
of the plateau, as in Fig. [5^. The depinning transition line F c shows a jump Av in 
the average v cm and a clear hysteretic behavior [16] as F crosses F c , for not too large 
values of K < K c . As can be expected, Av decreases to as K increases. Around 
K = K c the depinning transition is continuous. The specific point K = K c represents 
the zero-F crossing of the depinning line, which extends to negative F above K c . The 
precise value of K c depends on parameters such as t> ext and 7. As K approaches K c 
from above (no external force), v cm approaches ^plateau in a critical manner, as suggested 
in Fig. |5fa). This is detailed in Fig. [5](c,d), where the critical behavior is shown to be 
Av oc (K — Kc) 1 / 2 , the value typical of intermittencies of type I [14J. For K > K c , 
in fact the chain spends most of its time moving at v cm (t) ~ f plateau, except for short 
bursts at regular time-intervals r, where the system as a whole jumps ahead by a , 
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Figure 6. The average dissipated power per particle for the GM (r+,r_) = (</>, <fi 2 ) 
model as a function of the spring stiffness K. All dynamical parameters are as in 
Fig. [21 Two dotted lines indicate the plateau range. A dashed horizontal line marks 
the theoretical lowest dissipation limit of a rigid lubricant in a purely translational 
motion at speed v p i at cau = (1 - v cxt - 

i.e. an extra chain lattice spacing (collective slip). The characteristic time r between 
successive collective slips diverges as r a (K — K()~ x l 2 for K — > K c , consistent with 
the critical behavior of w. We verified that the u>-plateaus for more general values of 
r + and of r_ show the same kind of infinite stiffness, and a critical decrease of F c near 
the plateau edge, similar to that of Fig. [5h for the GM. 

We have investigated the tribological properties of the model in the plateau regime. 
As the dissipative term makes no work on particles that move at velocity ii = v ext /2, 
it is particularly convenient to evaluate the dissipated power in the "symmetric" frame 
of reference where the upper chain moves at velocity and the lower chain moves at 
velocity — -^f 1 - In this frame of reference, the instantaneous total power dissipated by 
the 7 term amounts to 

P = 5>;-( 2 7^) = 2 7 5> 4 2 = -^£ kin , (8) 

i i m 

where E^ n is the total kinetic energy of the particles in the symmetric frame of reference. 
This power fluctuates, but eventually, in the dynamical stationary state its average value 
(P) equals on average the mean power necessary to maintain the motion of the sliders 
[12] . Figure [6] shows the dissipated power per particle in a wide range of K values for 
the (4>, 4> 2 ) model. It is apparent that the friction decreases as the springs become stiffer 
and stiffer, where fluctuations are of decreasing amplitude. The plateau region does not 
appear to be especially conspicuous in this frictional decrease, except for a net drop at 
the end of the plateau, due to the lubricant moving closer and closer to the symmetric 
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velocity v ext /2. 
3. Conclusions 

The phenomena just described for a model ID system might have 2D counterparts, 
potentially observable in real systems. Nested carbon nanotubes [XT] , or confined one- 
dimensional nanomechanical systems [18] . are one of several possible arenas for the 
phenomena described in this paper. Though speculative at this stage, an obvious 
question is what aspects of the phenomenology just described might survive in 2D, 
where tribological realizations, such as the sliding of two hard crystalline faces with, 
e.g., an interposed graphite flake, are conceivable. Our results suggests that the lattice 
of discommensurations - a Moire pattern- formed by the flake on a substrate, could be 
dragged by the other sliding crystal face, in such a manner that the speed of the flake as 
a whole would be smaller, and quantized. This would amount to the slider "ironing" the 
solitons onward. Dienwiebel et al. [19J demonstrated how incommensurability may lead 
to virtually friction-free sliding in such a case, but no measure was obtained for the flake 
relative sliding velocity. Unlike our model, real substrates are not rigid nor ideal, but 
subject to thermal expansion and characterized by defects. Nevertheless the ubiquity of 
plateaus shown in Fig. [2J and their topological origin, suggests that these effects would 
not remove the phenomenon. A real-life situation with a distribution of differently 
oriented crystalline micro-grains, each possessing a different incommensurability, is also 
potentially interesting; each grain, we expect, will tend to stabilize a certain average 
CM velocity depending on its incommensurability. Other realizations or applications 
inspired by the physics described by our model might be accessible, notably in grain- 
boundary motion, in the sliding of optical lattices [20J or of charge- density- wave systems 

m- 
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